Scaling patterns of the suppression of 7r° yields in Au+Au collisions at 
$nn = 200 GeV: links to the transport properties of the QGP 
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Suppression measurements for neutral pions (it ) are used to investigate the predicted path length 
(L) and transverse momentum (pr) dependent jet quenching patterns of the hot QCD medium pro- 
duced in Au+Au collisions at ^/sjvjv = 200 GeV. The observed scaling patterns show the predicted 
trends for jet-medium interactions dominated by radiative energy loss. They also allow simple esti- 
mates of the transport coefficient q and the ratio of viscosity to entropy density n/s. These estimates 
indicate that the short mean free path (A) in the QCD medium leading to hydrodynamic-like flow 
with a small value of rj/s, is also responsible for the strong suppression observed. 



One of the important discoveries at the Relativis- 
tic Heavy Ion Collider (RHIC) at Brookhaven National 
Laboratory (BNL), has been the observation that high- 
Pt hadron yields are suppressed in central and mid- 
central A+A collisions when compared to the binary- 
scaled yields from p+p collisions This observation 
has been attributed to jet-quenching Q - the process by 
which hard scattered partons interact and loose energy 
in the hot and dense quark gluon plasma (QGP) pro- 
duced in the collisions. Subsequent to such interactions, 
the partons which do emerge, then fragment into topo- 
logically aligned hadrons (jets) which provide the basis 
for the 7T° suppression measurements. 

There is considerable current interest in the use of jet 
quenching as a quantitative tomographic probe of the 
QGP. Recent theoretical efforts have centered on inves- 
tigations of the energy loss mechanism for scattered par- 
tons which propagate through this medium. Two such 
mechanisms are; (i) scatterings off thermal partons in 
binary elastic collisions and (ii) Gluon bremsstrahlung 
with the Landau-Pomeranchuk-Migdal (LPM) 0] effect. 
The latter has been investigated via different formalisms 
H H B 0, H B ■ Studies of the relative importance 
of both jet quenching mechanisms have also been made 
[U EH 0, EH • To date, a conclusive mechanistic picture 
has not yet emerged. 

Initial quantitative studies with models which incorpo- 
rate the time evolution of the QGP medium via relativis- 
tic ideal (3+l)-dimensional hydrodynamical simulations, 
are currently underway [HI, [r|. However, the value of 
such studies rests heavily on accurate knowledge of the 
dominant mechanism/s for jet quenching. Therefore, it 
is important to pursue validation tests which can lend 
insight or provide a clear distinction between different 
energy loss mechanisms. 

The experimental probe commonly exploited for jet- 
quenching studies in AA collisions is the nuclear modifi- 
cation factor (Raa)', 
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FIG. 1: (Color online) Nuclear modification factor (Raa) for 
7r°s reproduced from Ref. [TJ- Error bars indicate statisti- 
cal and px-uncorrelated errors; boxes indicate px-correlated 
errors. The single box around -Raa=1 on the left is the er- 
ror due to JVcon; the single box on the right is the overall 
normalization error of the p+p reference spectrum. 



collisions and (Taa) is the nuclear thickness function av- 
eraged over the impact parameter range associated with 
a given centrality selection 



(Taa) = 



J T AA (b) dh 



/(I 
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where er™ is the particle production cross section in p+p 



The corresponding average number of nucleon-nucleon 
collisions, (N co u) = <j™ el ' (T aa) , is routinely obtained via 
a Monte-Carlo Glauber-based model calculation [l8l. Il9j|. 

In this letter we use Raa measurements to perform 
validation tests which addresses the question of whether 
or not medium induced gluon radiation is a dominant 
mechanism for jet-energy loss. 
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The measurements employed for these tests are the 
recently published Raa{pt) and Raa{^-4>tPt) data for 
7T° [l7L I20I]. Here, A</> is the azimuthal angle relative to 
the reaction plane. A subset of these data is shown as a 
function of px for several centrality selections in Fig. [1] 
The minimum bias results in the top left panel indicate 
that, for 2.5 < pt < 5 GeV/c, suppression increases with 
Pt to the value -Raa ~0.3. By contrast, a much weaker 
Pt dependence is observed for pt > 5 GeV/c, with a hint 
that the suppression decreases as pr increases. The same 
trends are evident for all centrality selections spanning 
central and mid-central collisions, albeit with different 
absolute magnitudes. These trends provide an important 
constraint for our tests, as discussed below. 

To perform validation tests on the data for medium 
induced gluon radiation, we use the "pocket formula" of 
Dokshitzcr and Kharzeev (DK) Q. This formula gives 
the quenching of the transverse momentum spectrum for 
jets produced from scattered light partons [8|] as; 



Raa(pt,L) ~ exp 
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where a s is the strong interaction coupling strength, Cp 
is the color factor, L is the path length [of the medium] 
that the parton traverses and q is the transport coefficient 
which reflects the squared average transverse momentum 
exchange per unit path length, between the medium and 
the parton. 

Here, there are two essential points. First, Eq. [1] is 
derived with an explicit assumption that the mechanism 
for energy loss is medium induced gluon radiation [2l| . 
Second, this equation gives specific testable predictions 
for the dependence of Raa{pt, L) on L and pt for pt > 
5 GcV/c. That is, In [Raa(pt, L)] should scale as L and 
1/a/(pt) respectively (ie. In [Raa(pt, L)] should show 
a linear dependence on both L and 1/ \/~(pt) if medium 
induced gluon radiation is indeed the dominant energy 
loss mechanism). 

To study the L dependence, the transverse size of the 
system R was used as an estimate for the angle averaged 
path length L, as well as to determine its value in- (L x ) 
and out- (L y ) of the reaction plane. For each central- 
ity selection, the number of participant nucleons 7V part , 
was estimated via a Monte-Carlo Glauber-based model 
[HI Eil ■ The corresponding transverse size R was then 
determined from the distribution of these nucleons in the 
transverse (a;, y) plane via the same Monte-Carlo Glauber 
model: 




where a x and a y are the respective root-mean-squarc 
widths of the density distributions; here, averaging is per- 
formed over configurations. For these calculations, the 
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2: (Color online) In [Raa(pt, L)) vs. l/y^pr) (pt > 
5 GeV/c) for centrality selections of 0-5% (a) and 20-30% 
(b). Error bars are statistical only. pT-correlated systematic 
errors are about 12% and the systematic errors associated 
with JVcoii and the overall normalization of the p+p reference 
spectrum are ~ 7% and 10% respectively. The dot-dashed 
curve in each panel is a fit to the data (see text). 



initial entropy profile in the transverse plane was assumed 
to be proportional to a linear combination of the num- 
ber density of participants and binary collisions [13, [HI . 
The latter assures that the entropy density weighting is 
constrained by hadron multiplicity measurements. 

The results from our validation tests for medium- 
induced gluon radiation are summarized in Figs. [5] - 3] 
Figs. [2] (a) and (b) show plots of In [Raa(pt, L)] vs. 
l/y/lpr) for pt > 5 GeV/c and centrality selections of 
0-5% (L=1.90 fm) and 20-30% (L=1.45 fm) respectively. 
The dot-dashed curves in these figures represent a lin- 
ear fit to the data. They show that, although there is 
some scatter, the trend of the data is compatible with 
the l/y/lpx) dependence predicted by Eq. [T] ft is worth 
mentioning here that a similar 1/\/(pt) dependence is 
observed for the range 2.5 < pt < 5 GeV/c, but with a 
different (positive) slope. 

Figure[3]shows the L dependence of In [Raa(Pt, L)\ for 
two different pt selections as indicated. The dot-dashed 
curve shows a linear fit to the data for pt > 5 GeV/c and 
L > 0.65 fm, ie. the two data points corresponding to the 
two most peripheral collision ccntralities were excluded 
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FIG. 3: (Color online) In [Raa(pt, L)] vs. L for two px se- 
lections as indicated. Error bars are statistical only. The 
systematic error associated with -/V co n ranges from ~ 7% in 
central collisions to ~ 30% in peripheral collisions; the sys- 
tematic error resulting from overall normalization of the p+p 
reference spectrum is about 10% (cf. FigfTJ. The dot-dashed 
curve is a linear fit to the data set for pr > 5 GeV/c (see 
text). 



from the fit (note the large N co \\ systematic errors for 
peripheral collisions). Fig. [3] shows that the data trends 
for both pt selections (px > 5 and px > 10 GeV/c), 
validates the linear dependence on L predicted by Eq. [T] 
The efficacy of this scaling pattern is not significantly 
influenced by systematic errors. 

Further evidence for this linear dependence is also ev- 
ident in Fig. [4] where we show results for Raa measure- 
ments in- and out of the reaction plane for 7r°s with 
5 < p T < 8 GeV/c [2(|. Panel (a) shows that, for the 
same number of participants. 7r° suppression is larger out- 
of-plane than it is in-plane. However, panel (b) shows 
that, when plotted as \n[RxA{PT,L x ^ y )] vs. L x and L y , 
the in-plane and out-of-plane data show a single linear 
dependence on i x y . This same dependence is found for 
all other reaction plane orientations for pt > 5 GeV/c. 
This confirms that the azimuthal anisotropy of particle 
yields (at higher) stem from the path-length dependence 
of jet quenching. The dashed curve (fit to the data) in 
the figure reinforces the predicted linear dependence of 
these data on the path length (cf. Eq. [T|) . 

The fits to the data in Fig. [3] (for the two px cuts) indi- 
cate an intercept L w 0.6 fm, which suggests a minimum 
path length requirement for the initiation of ir° suppres- 
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FIG. 4: (Color online) (a) Raa(pt, N palt ) vs. jV part for 7r°'s 
detected in- and out of the reaction plane. The data are taken 
from Ref. [13] for the selection 5 < pr < 8 GeV/c. Error 
bars are statistical only, (b) \ti[Raa(pt, L)] vs. I/ X , y , the 
path length in- and out of the reaction plane. The systematic 
errors for both (a) and (b) are similar to those of Fig. [3] The 
dashed curve is a linear fit to the data. 



sion. Such a requirement is akin to the plasma formation 
or cooling times proposed in Refs. 0, H^l- 

The corresponding slopes (for px > 5 and px > 10 
GeV/c) indicate similar magnitudes with a hint of a small 
decrease in slope with increasing (pt)- This relatively 
weak px dependence (which is also reflected in Fig. [2]), 
could be an indication that q has a dependence on px, eg. 

q = q~Q fp-J where A is a small fractional power (26| . 

To obtain an estimate of the magnitude of q, we use 
Eq. [T]in conjunction with the slope, —1.26 ± 0.06 fm -1 , 
extracted from Fig. for px > 10 GeV/c. This gives 
the value q « 0.75 GeV 2 /fm for the values a s = 0.3 
0, C F = 9/4 @, S3] and C = n = 8.1 ± 0.05 0. 
This estimate of q, which can be interpreted as a space- 
time average, is comparable to the recent estimates of 
~ 1— 2 GeV 2 /fm obtained from fits to hadron suppression 
data (for the most central Au+Au collisions) within the 
framework of the higher twist (HT) expansion [28l. Sil 
and the Gyulassy-Levai-Vitev (GLV) scheme [30l. l3l|. It 
is, however, much less than the value recently extracted 
via the approach of Arnold, Moore and Yaffe (AMY) 
[ill SH and that of Armesto Salgado and Wiedemann 

(ASW) plum. 

The ratio of the shear viscosity (rj) to the entropy den- 
sity (s) can also be estimated as |34| ; 

n T 3 
- w 1.25 — 

s q 

where T is the temperature. This estimate can be com- 
pared to the value of Air^ m 1.3 ± 0.3 extracted from 
flow data SH for T - 220 MeV H3|. It is noteworthy 
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that if we use this same temperature in conjunction with 
our extracted value for q (in the above equation) , we ob- 
tain a strikingly similar estimate for rj/s. We conclude 
therefore, that the relatively short mean free path in the 
plasma |23| which drives hydrodynamic-like flow with 
small shear viscosity is also responsible for the strong 
jet quenching observed. 

In summary we have performed validation tests of the 
scaling properties of jet quenching to investigate the dom- 
inant mechanism for jet-energy loss. Our tests confirm 
the weak 1/\/(j>t) dependence, as well as the linear de- 
pendence on path length (L and L x ^ y ) predicted by Dok- 
shitzer and Kharzeev, for jet suppression dominated by 
the mechanism of medium-induced gluon radiation in a 
hot and dense QGP. The quenching patterns indicate a 
minimum path length requirement for the initiation of 
7T° suppression (corona), suggests a possible pt depen- 
dence for q and gives the estimate q ~ 1 GeV 2 /fm which 
is comparable to that obtained within the framework of 



the HT and GLV models. This estimate also indicates 
a small value of rj/s which is in good agreement with 
that obtained via flow measurements. The present study 
will have to be extended to the heavy quark systems to 
arrive at an even more definitive conclusion about the 
role of a perturbative radiative energy loss mechanism. 
Nonetheless, these results will undoubtedly provide im- 
portant model constraints in future attempts to use jet- 
quenching as a quantitative tomographic probe of the 
QGP. 



Acknowledgments 

We thank D.E. Kharzeev for crucial and interesting 
discussions. This work was supported by the US DOE 
under contract DE-FG02-87ER40331.A008. 



[1] K. Adcox et al. (PHENIX), Phys. Rev. Lett. 88, 022301 

(2002), nucl-ex/0109003. 
[2] M. Gyulassy and X.-n. Wang, Nucl. Phys. B420, 583 

(1994), nucl-th/9306003. 
[3] A. B. Migdal, Phys. Rev. 103, 1811 (1956). 
[4] R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne, 

and D. Schiff, Nucl. Phys. B483, 291 (1997), hep- 

ph/9607355. 

[5] A. Kovner and U. A. Wiedemann (2003), hep- 
ph/0304151. 

[6] B. G. Zakharov, JETP Lett. 63, 952 (1996), hep- 
ph/9607440. 

[7] M. Gyulassy, P. Levai, and I. Vitev, Nucl. Phys. B594, 

371 (2001), nucl-th/0006010. 
[8] Y. L. Dokshitzer and D. E. Kharzeev, Phys. Lett. B519, 

199 (2001), hep-ph/0106202. 
[9] X.-N. Wang and X.-f. Guo, Nucl. Phys. A696, 788 

(2001), hep-ph/0102230. 
[10] A. Majumder and B. Muller, Phys. Rev. C77, 054903 

(2008), 0705.1147. 
[11] M. G. Mustafa and M. H. Thoma, Acta Phys. Hung. 

A22, 93 (2005), hep-ph/0311168. 
[12] A. Adil, M. Gyulassy, W. A. Horowitz, and S. Wicks, 

Phys. Rev. C75, 044906 (2007), nucl-th/0606010. 
[13] B. G. Zakharov, JETP Lett. 86, 444 (2007), 0708.0816. 
[14] T. Renk, Phys. Rev. C76, 064905 (2007), 0708.4319. 
[15] G.-Y. Qin et al., Phys. Rev. Lett. 100, 072301 (2008), 

0710.0605. 

[16] S. A. Bass et al., Phys. Rev. C79, 024901 (2009), 
0808.0908. 

[17] A. Adare et al. (PHENIX), Phys. Rev. Lett. 101, 232301 
(2008), 0801.4020. 

[18] M. L. Miller, K. Reygers, S. J. Sanders, and P. Stein- 
berg, Ann. Rev. Nucl. Part. Sci. 57, 205 (2007), nucl- 



ex/0701025. 

[19] B. Alver et al., Phys. Rev. Lett. 98, 242302 (2007). 
[20] S. S. Adler et al. (PHENIX), Phys. Rev. C76, 034904 

(2007), nucl-ex/0611007. 
[21] Radiative energy loss is expected to grow quadratically 

with path length while elastic energy loss is expected to 

grow linearly (see for example Ref. |14[h 
[22] T. Hirano and Y. Nara (2009), 0904.4080. 
[23] R. A. Lacey, A. Taranenko, and R. Wei (2009), 

0905.4368. 

[24] V. S. Pantuev, JETP Lett. 85, 104 (2007), hep- 
ph/0506095. 

[25] J. Liao and E. V. Shuryak, Nucl. Phys. A775, 224 (2006), 

hep-ph/0508035. 
[26] J. Casalderrey-Solana and X.-N. Wang, Phys. Rev. C77, 

024902 (2008), 0705.1352. 
[27] Here, it is assumed that ir°s are dominantly produced 

via the fragmentation of gluons. 
[28] H. Zhang, J. F. Owens, E. Wang, and X.-N. Wang, Phys. 

Rev. Lett. 98, 212301 (2007), nucl-th/0701045. 
[29] A. Majumder, C. Nonaka, and S. A. Bass, Phys. Rev. 

C76, 041902 (2007), nucl-th/0703019. 
[30] M. Gyulassy, P. Levai, and I. Vitev, Phys. Rev. Lett. 85, 

5535 (2000), micl-th/0005032. 
[31] A. Majumder, J. Phys. G34, S377 (2007), nucl- 

th/0702066. 

[32] G.-Y. Qin et al., Phys. Rev. C76, 064907 (2007), 
0705.2575. 

[33] U. A. Wiedemann, Nucl. Phys. B582, 409 (2000), hep- 
ph/0003021. 

[34] A. Majumder, B. Muller, and X.-N. Wang, Phys. Rev. 

Lett. 99, 192301 (2007), hep-ph/0703082. 
[35] A. Adare et al. (PHENIX) (2008), 0804.4168. 



